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We study the effect of the voltage bias on the ferromagnetic phase transition in a one-dimensional 
itinerant electron system. The applied voltage drives the system into a nonequilibrium steady 
state with a non-zero electric current. The bias changes the universality class of the second order 
ferromagnetic transition. While the equilibrium transition belongs to the universality class of the 
uniaxial ferroelectric, we find the mean-field behavior near the nonequilibrium critical point. 



When external fields drive a many-body system into a nonequilibrium steady state, a parameter change may result 
in a transition between two nonequilibrium phases. Such nonequilibrium transitions were observed in 2D electron 
gases driven with microwave radiation 0,0 , granular matter and many other systems. Scaling and universality 
hold near second order phase transitions both in and far from thermal equilibrium jj]. Hence, one can describe them 
with rcnormalization group. This helped to achieve significant progress in the theory of nonequilibrium classical phase 
transitions Q. 

The possibility to achieve a high degree of isolation from environment for ultracold atoms [5| opens a way to 
study quantum dynamics in many-body systems far from equilibrium. This has stimulated interest to nonequilibrium 
quantum phase transitions of strongly interacting bosons [6(. Another type of nonequilibrium quantum transitions 
is possible in itinerant electron systems driven into states with nonzero electric current. In this Letter, we address a 
nonequilibrium ferromagnetic transition of ID itinerant fermions. Since this transition exists only at zero temperature, 
it cannot be considered as a classical phenomenon, in contrast to the decay of supercurrents of interacting bosons 
that can be understood as a classical modulational instability p. 

Equilibrium magnetic phase transitions of itinerant electrons have attracted much attention 0> 0, 0> E3- The 
quantum critical behavior at these transitions is a subtle and challenging problem due to the coexistence of two types 
of soft modes. In addition to the order parameter fluctuations, the theory must account for gajdess fermions which 
induce effective long-range interactions between the order parameter fields in different points Most studies have 
focused on itinerant systems in two and three dimensions. According to the Lieb and Mattis theorem there is 
no ferromagnetism in a wide class of one-dimensional models with spin-independent interactions. Yet ferromagnetic 
quantum phase transitions were found in many one-dimensional itinerant systems both with and without magnetic 
anisotropy. In particular, recent numerical simulations [l2T | discovered a quantum critical point in a tight-binding 
model with next nearest neighbor hoping. It was also argued that the 0.7 anomaly in transport properties 
of quantum wires flU , Il4j could be explained by spontaneous magnetization of itinerant electrons. Recently, the 
equilibrium quantum phase transition of itinerant electrons in ID was theoretically investigated in Refs. Il5llla . 

In this Letter we study the effect of a voltage bias on the quantum ferromagnetic transition in a one-dimensional 
itinerant electron system. A low voltage cannot significantly affect the magnetization far from the transition. However, 
it is a relevant perturbation in the critical domain and changes the universality class of the transition. While the 
equilibrium second order transition 0] belongs to the universality class of the uniaxial ferroelectric 0, llil H^ |. 
we find that the nonequilibrium critical behavior is mean-field like, i.e. the magnetization m ~ y^r, where r is the 
distance from the critical point. 

The paper is organized as follows: first we introduce the model, then briefly discuss the phase transition within the 
mean-field Stoner approach, and finally show that the mean-field approximation gives the exact critical behavior. 
We consider the Tomonaga-Luttinger model with the Hamiltonian 



H = J dxj Y,bPln(xX-Mdx)i>Mx) + Tpt L (x)e(-ihd x )i> aL (x)} + XJ / dyK ° 



'(x-y)pa(x)pa'(y)>, (1) 



where ip]y R and ip a L are the creation operators for right- and left-moving electrons with the spin z-component a = 
±1/2, = V^fl + Vvz gives the conventional electron creation operator, p a — ip^tpu is the electron density, and 
Kaa' (x — y) the interaction strength. We assume that the long-range Coulomb interaction is screened by the gates 
so that K aa >(x — y) decreases rapidly for large (x — y). The spin dependence of K aa i determines the symmetry of 
the model. As in Ref. 0we will focus on the technically simplest case of the Ising symmetry Z2 x U(l) (U(l) is 
responsible for the conservation of S z ] Z% describes the symmetry with respect to the a — > — a transformation). Note 
that the nonlinearity of the spectrum e(k) is relevant near the critical point. As discussed in Ref. ^^this nonlinearity 
is important for the very existence of the quantum critical point. 
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The equilibrium transition in the model was considered in Ref. We will study what changes in the presence 
of a voltage bias. We will show that the non-equilibrium problem is equivalent to a certain equilibrium model which 
lacks the time-reversal symmetry. The absence of this symmetry is related to the fact that the "Fermi velocities" of 
the right- and left-moving electrons are different. The critical behavior of the latter model can be found exactly. We 
will show that the diagrams for the magnetization are free from infra-red divergencies. Hence, the system belongs to 
the mean-field universality class. 

Before we develop a systematic approach let us briefly discuss the simplified Stoner picture. Let the number of 
the right/left-movers with spin a be No + Nr/lc where No denotes the number of the states with the negative 
energy e(k). The Stoner model is described by the Hamiltonian H = J2 a J2d=r. l U(Noa) — KM 2 , where U(n) = 

I^N e (^)|f ' an d the magnetization M — [N R [i/ 2 ] + N L [i/ 2 ] — ^Vr[-i/2] — ^Yl[-i/2])/2- The numbers of the 
right- and left-movers are integrals of motion. The chemical potentials of the right- and left-movers (J.r/l — ±eV/2 
are different. Hence, at zero temperature the system is in the ground state of the effective Hamiltonian H' = 
H - p R N R - p, L N L = H - V(N R [ X /2] + Nr[_i/ 2 ] - N L[1 / 2 ] - N L ^ 1/2] )/2. From the energy minimum condition 
dH 1 / dNpu = 0we find that the maximal energies e(fc) of the right-moving particles with the opposite spins differ by 
2AE = K f_£ E dE[p(E + V/2) + p(E - V/2)], where p(E) is the density of states. Expanding p(E ± V/2) in powers 
of E, V one finds AE[1 - 2K(p(0) + p"(0)V 2 /8)} - KAE 3 p"(0)/3 = 0. This equation describes a second order phase 
transition, if p"(0) < 0. The voltage bias shifts the transition to the greater interaction strength K = K c , AK C ~ V 2 . 
The magnetization exhibits the mean-field behavior M ~ AE ~ \[K — K c . We will see below that the mean field 
prediction for the critical behavior of the magnetization is valid beyond the Stoner approximation. 

In the Stoner picture, electrons move in a spin-dependent self-consistent field without backscattering. Hence, the 
electric current is I = 2e 2 V/h. The spin current I s = 0. 

In a systematic approach^ we need to include the leads into the model to describe the voltage bias. We use the 
standard model |2Jj, |2jJ, |22j, |2jj for the Fermi-liquid leads adiabatically connected to the wire. We assume that the 
Hamiltonian Q is applicable for \x\ < L only, where 2L is the wire length. At large \x\ the interaction strength 
K(x — y) is zero. We assume that at the initial moment of time t = — oo the interaction K aa * = 0. Thus, the numbers 
of the left- and right-movers are integrals of motion at t = — oo. Hence, the initial state can be described in terms 
of two different chemical potentials [ir — eV/2 and pl — ~eV/2 for the right- and left-moving electrons, where V 
is the voltage bias. At later times t > — oo the interaction is gradually turned on. It is convenient to switch to 
the interaction representation so that the Hamiltonian becomes H' = H — plNl — PrNr, where Nl and Nr are 
the numbers of the left- and right-movers. The initial state is the ground state of H' at t = — oo. The interaction 
representation introduces time dependence into the electron creation and annihilation operators. Hence, the evolution 
of the system should be described with the nonequilibrium Keldysh formalism [2ll I22I Eij . Eif . The bare Keldysh 
Green functions are determined by the ground state of the Hamiltonian H' at t = —00. We will see however, that the 
time-dcpcndcncc drops out from the problem and a simpler equilibrium technique can be used. 

The bosonization approach |25| allows one to reduce one-dimensional problems with strong electron interaction to 
models of weakly interacting bosons which can be studied perturbatively in the interaction strength. Following the 
standard bosonization procedure we define four Bose-fields 4> a R and such that 

ip aR — exp(i0 CTfl + ikfRX - Wt/2);ipaL ~ exp(-z0 CTL - ik fL x + iVt/2), (2) 

where kfR and kfL play the role of the effective "Fermi momenta" and are proportional to the average electron 
density of the right- and left-movers, p L/R = Y^aP^L/R = d x<PaL/R + 2fc /L/jR )/(27r), (pl/r) = k fL / H /ir. We 
omit Klein factors in Eqs. (0 since they are not important for the following calculations. Two wave vectors ^/l/_r 
differ since the applied voltage bias leads to the different densities of the right- and left-moving electrons. In most 
cases this difference can be neglected since the voltage is much smaller than the bandwidth. However, it is important 
in our problem. The bosonized action contains two types of terms. The first type is a product of derivatives of the 
Bose-fields with respect to t and x ("forward scattering" terms). The second type includes "backscattering" terms 
proportional to exp(zn[(/> c /2 + (kfR + kfi,)x — Vt] + il<p s /2), where the charge and spin fields <f> c = ^2 a {4>aR. + 4>ol) 
and 4> s — J2 a 2cr((f> a -R + 4> a L) are related to the charge and spin densities as p — e(d x <f) c + 2kfL + 2fc/i?)/(27r) and 
m = d x (f> s / '(47r); n and I are integers, (n + l) is even. The dual fields 9 C = J2 a {4>a-R — 4><tl) and 6 S = J2 a 2a(ip a -R — 4> a L) 
do not appear in the exponents due to the conservation of the charge and the z-component of the magnetization. The 
terms of the second type with n ^ are suppressed by their rapidly oscillating coordinate dependence. Hence, as 
usual in bosonization they are not important for the low-energy physics [2f| which determines the critical behavior. 
We will also see that the terms with n — 0, 1 7^ are irrelevant near the phase transition. Thus, only the terms of 
the first type should be taken into account. This means that the terms which should be included in the action do not 
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depend on time explicitly. Hence, it becomes possible to use the equilibrium imaginary time formalism instead of the 
non-equilibrium Keldysh technique |2]J, |22j, 1 2.1 |24| . 

The Euclidean action contains four interacting chiral fields 4> a R/L'- 

L= dxdrl^ , ^[d x (j> aR (id T +v R d x )(f) rTR +d x (j> a L(-idT+VLd x )(t) a L}+ u aa 'DD'd x <i>aDd x <l><T>D>+- ■ ■ 

cr a.a' = ±l/2;D.D'=R,L 

(3) 

where dots denote higher order terms. Note that the "Fermi velocities" Vr/l — de ( k ~^t R / L '> are different for the right - 
and left-movers. We rewrite the action in terms of </) c ,(f) s ,9 c ,9 3 and integrate out 8 S and 8 C . The symmetry with 
respect to the magnetization reversal results in the spin-charge separation, i.e. the quadratic part of the resulting 
action is the sum of a contribution which depends on <j) c only and a contribution that depends on <j) s only: 



L^=L^(4> S ) + L^{4> C ). (4) 

In general, the spin-charge separation does not extend beyond the terms quadratic in the Bose-fields. However, similar 
to Ref. the interaction between the spin and charge degrees of freedom is not important near the phase transition. 
One can show that this interaction does not change the critical behavior. Thus, we can focus on the case when 
Eq. Q is valid in all orders and hence the magnetization depends only on the spin part of the action L s ((j) s ). A 
straightforward modification of our approach demonstrates that our results are independent of this assumption (cf. 
Ref. El). 

After an appropriate rescaling of the time variable, the action takes the following form: 



L = J ' dTdx^[(d T s ) 2 + 2icd T <j) s d x cj )s +r{d x ^ s f] + ...j, (5) 

where the dots denote terms with higher powers of <f> s and/or gradients. The terms with an odd (even) number of 
time derivatives are imaginary (real). The action JSJ) differs from the standard expression 25] due to the presence of 
the second imaginary term. Technically this term appears because two "Fermi velocities" v R and vl are different. 
Its presence reflects the fact that the action (0 lacks time-reversal symmetry. In most problems the second term can 
be neglected. Indeed, the constant c is proportional to the voltage bias which is much smaller than the bandwidth. 
In particular, the second term is not important far from the transition since it is smaller than the sum of the first 
and third terms in (0). However, at the critical point, r vanishes and hence the second term is not negligible. At 
negative r the quadratic part of the action JSJ) becomes unstable. This signals a phase transition. To make sure that 
the system is stable at negative r one needs to take into account the contribution of the form £2,2 = b(d 2 (f) s ) 2 and 
nonlinear terms such as Lx,4 = u (d x (f>s}^ denoted by dots in Eq. JSJ). L 2i 2 emerges from the coordinate dependence 
of the interaction K{x — y), Eq. Q, [lj|. L 1;4 describes the nonlinearity of the spectrum e(fc) |26j |. 

In the critical point the quadratic part of the action is thus proportional to {d T <f) s ) 2 + 2icd T (f) s d x 4> s + b{d 2 (j) s ) 2 . 
At each length scale 1/k the poles of the propagator determine two frequency scales U!% ~ k and U2 ~ k 3 . Hence, 
similar to other models of itinerant electron systems 0, the choice of scaling dimensions is subtle. Indeed, if the 
scaling dimension of the coordinate [x] = — 1 then two different frequency scales correspond to [t] = — 1 and [t] = —3. 
However, for any of these two choices, power counting yields a negative or zero scaling dimension for all possible 
non-linear operators. This suggests the mean-field behavior. Below we show that the critical behavior is indeed 
mean- field-like without logarithmic corrections. 

The idea of our calculations can be understood from the following example. According to the Ginzburg criterion, 
deviations from the mean-field theory emerge if the diagram depicted in Fig. 1 diverges in the infra-red limit, i.e. when 
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incoming momenta and frequencies are low. The diagram reduces to the integral I = J dkdujk 4 /[uj 2 + 2icuik + bk 4 ] 2 . 
We substitute uj = ui' — ick and shift the integration contour over uj' onto the real axis. This gives I = J dkduj' k 4 /[uj' 2 + 
c 2 k 2 + bk 4 ] 2 . The latter integral converges. 

To calculate the magnetization we add to the action a weak magnetic field term, —h q iqcj) s (q)/(4:'K), where 4> s (q) 
denotes a Fourier component of the spin Bose-field. The uniform field corresponds to the limit q —> 0. At ft, = 
the magnetization is a non-analytic function of r but for any finite magnetic field the magnetization m depends on r 
analytically [2?1 Eij . Thus, the critical behavior of the magnetization at small negative r can be determined with the 
following three steps [13] : 1) calculate m = (d x (f> s ) / (An) at r > using the perturbation theory about the cj) s = 
state; 2) analytically continue m(r) to negative r; 3) take the limit h — > 0. At r > we get the standard result |27j 

f h °° 1 1 

(d x <t> s {x = 0)> = Cj— - (2n-iy T2n ^ Ql = ■ ' • = 92,1 = w i = • ' • = w 2« = 0){8*M* = O)) 2 "- 1 j; (6) 



C= lim limg 2 G= l/r;G= \/[auj 2 + 2i£ojq + fq 2 ], (7) 

where T2 n are the proper vertices at zero momentum and frequency, and a, c and f denote the renormalized coefficients 
of the Green function G. The formulas for the proper vertices contain integrals of the products of the expressions 
of the form pip2G(uj,q), where pi,P2 = q. It is convenient to change variables: uj' = ui + icq/a. Then G — ► G' — 
\j\aui' 2 + (c 2 /a + r)q 2 ). At f > the integration path for integrals over dui' can be deformed into the real axis. 
(One can show this, e.g., using the representation G' = J °° ds exp{— s[auj' 2 + (c 2 /a + f)q 2 ]} in Feynman integrals.) 
The resulting integrals have no infrared divergences at f = 0. Expanding G' in powers of f one can see that those 
integrals are analytic functions of f at small f. Thus, all coefficients of Eq. are finite at f = and can be 
analytically continued to negative f. Eq. JBJ has the structure of the mean-field equation with the free energy 
F = —lira + f(47rm) 2 /2 + ^2 ^2n{r)(^TTm) 2n / (2n)\. Depending on the coefficients it may describe both first and 
second order phase transitions. A second order transition is possible for positive T±{f = 0). In this case the average 
magnetization follows the mean-field laws: 



m{f <0,h = Q)~V-r; dm ( h = > f >°) „ l Mr = 0, h) ~ (8) 

dh r 

The magnetization correlation function does not have the standard mean-field form. This is not surprising since 
in Luttinger liquids the correlation length is infinite for any r and not only in the critical point j^. Using the 
perturbation theory about the state with d x (j) s = 47rm(r) we find the large-distance behavior of the correlation 
function 



i(d x (j> a {x = Q,t = Q)-^m){d x (j>Jx = X,t = Q)-^m)) = -d 2 x \ r- exp(ikX)G(uj, k) - 1/A 2 . (9) 

J Air 2 

The coefficient before 1 / X 2 is a non-analytic function of f near the phase transition. 

So far we did not verify that cos(Z</> s ) is an irrelevant operator. This can be easily done with the renormalization 
group. At each step we represent <f> s = Airmx + <fif + <fif as the sum of the average Airmx, slow fluctuating part <f>f, 
and fast part = f^du; f bA< ^ <A dk exp(iwT + ikx)4* s (u>,k), where A is the ultra-violet cutoff and b <C 1. We 

integrate out the fast modes <j>f and make the rescaling k —> k/b, uj — > u)/b. One can easily find the renormalization 
of the cosine term in the zero-loop order. We see that cos(l<fi s ) — > b~ 2 cos(47rZmx/6 + l(f>f)ex-p(—l 2 ((f>f 2 )/2), where 
(4>f 2 ) ~ J du> J bA< ^ k ^ <A dkG(uj, k) ~ —(In 6)/ \/c 2 + a.K,\r\. Thus, the scaling dimension of the operator cos(/0 s ), 

2 — const/ •v/c 2 + a«|f|, is negative for small c ~ V and f. Hence, the operator is irrelevant near the phase transition. 

The above perturbative argument applies only if the coefficients wi before cos(l<p s ) are small. This explains why 
our results are valid in the case of the Ising symmetry Z2 x U(l) but not for the Heisenberg symmetry group even 
though the bosonized action has the same structure for both symmetries. Indeed, in the Heisenberg case w\ and r 
are related and cannot be simultaneously small |25j. 

The equilibrium ferromagnetic transition is described by the action L = (d T cf) s ) 2 /2+r(d x (f> s ) 2 /2+b(d 2 c/) s ) 4 +u(d x 6 £ ) 4: 
[lij . After the substitution m = d x 4> s /{Ait) this action reduces to the free energy of the uniaxial ferroelectric [l7j . 
Interestingly, the same universality class emerges in a number of unrelated non-equilibrium classical problems with 
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conserved order parameter 18, |l9j and in equilibrium quantum problems [29|. Long range dipole forces reduce the 
critical dimension of the uniaxial ferroelectrics to 3. This is equivalent to 2 in the quantum problem. Non-equilibrium 
effects are known to suppress fluctuations in classical systems 0, |3(|- Our results show that the same tendency 
exists in quantum itinerant systems with conserved S z where the mean-field behavior is possible in ID. On the other 
hand, non-equilibrium systems whose order parameter does not conserve exhibit typically the same critical behavior 
as corresponding equilibrium systems [3l|. 

In conclusion, we have shown that the voltage bias modifies the critical behavior of one-dimensional itinerant 
electrons near the ferromagnetic transition. The exact critical exponents are given by the mean field theory. 
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